
    

    Results.theta = array(NA, c(10, 2, 5))
    Results.phi   = array(NA, c(10, 2, 6))
    Results.gop   = matrix(NA, 10, 2)
    Results.nll   = matrix(NA, 10)

    for(seed in c(1:10)){
        load(paste("./cluster/submitted092616/RDataFiles/092616.mle.seed",seed,".RData",sep=""))   
        
        Results.theta[seed,,] = mle$theta.coef
        Results.phi[seed,,]   = mle$phi.coef
        Results.gop[seed,]    = mle$gop.coef
        Results.nll[seed]     = mle$nll
    }
    
    Results.nll
    Results.theta
    
    # load("./cluster/submitted092116/RDataFiles/v1/092116.mleREP1.RData")
    # 
    # nLL("par.update", parn = "not", mle, data, eta.coef, cnt, rep=500)
    # 